The dynamics of colloids in a narrow channel driven by a non-uniform force 
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Using Brownian dynamics simulations, we investigate the dynamics of colloids confined in two- 
dimensional narrow channels driven by a non-uniform force Fdr(y)- We considered linear-gradient, 
parabolic and delta-like driving-force profiles. This driving force induces melting of the colloidal 
solid (i.e., shear-induced melting), and the colloidal motion experiences a transition from elastic to 
plastic regime with increasing F^r- For intermediate Fdr (i-e., in the transition region) the response 
of the system, i.e., the distribution of the velocities of the colloidal chains Vi{y), in general does 
not coincide with the profile of the driving force Fdr{y), and depends on the magnitude of Fdr, 
the width of the channel and the density of colloids. For example, we show that the onset of 
plasticity is first observed near the boundaries while the motion in the central region is elastic. This 
is explained by: i) (in)commensurability between the chains due to the larger density of colloids 
near the boundaries, and ii) the gradient in Fdr- Our study provides a deeper understanding of 
the dynamics of colloids in channels and could be accessed in experiments on colloids (or in dusty 
plasma) with, e.g., asymmetric channels or in the presence of a gradient potential field. 

PACS numbers: 82.70.Dd, 64.60.Cn, 64.70.Rh 



I. INTRODUCTION 

During the last decade there has been a growing in- 
terest in the research of physical properties of colloidal 
systems. This interest is by part due to perspectives of 
practical use of colloids, e.g., in biology, medicine, meteo- 
rology, food production, etc. On the other hand, colloids 
serve as a convenient model system for studying, e.g., 
phase transitions diffusion p?], or commensurate- 

incommensurate transitions [5,]. Typical dimensions of 
colloids are in the micrometer regime and their dynamics 
is governed by a time scale in the microsecond regime and 
therefore is suitable for direct observation in real space 
and time. Furthemore, it has been possible to tune the 
inter-particle interaction potential which opens thereby a 
wide area for research of fundamental properties of clas- 
sical systems. 

In Refs. [1, 0, [1] structural properties of mag- 
netic mono-colloidal mixtures confined in narrow two- 
dimensional (2D) channels were studied under equilib- 
rium conditions. Structural deviations from an infinite 
2D crystal and related to it oscillations in structural 
properties were predicted. Transport properties of para- 
magnetic mono-colloidal mixtures in narrow channels be- 
ing in non-equilibrium but under stationary and homo- 
geneous external conditions (in the gravity field) were 
investigated in Refs. . It was shown that gravita- 

tion action results in the occurrence of a density gradient 
along the channel, that in its turn leads to a gradient in 
the number of chains. The latter was predicted for a 
driven colloidal system in the presence of a constriction 
in the channel [lll|. Besides, authors of Refs. also 
studied the relation between velocities of colloidal mo- 
tion, diffusive behavior and the self-organized order in 
the system. 



The transition from a hexagonal to a chain-like ordered 
phase was studied in Ref. [13] for Yukawa particle s (appl i- 
cable for charged colloids and dusty plasma [Tsl. [T4.]l5j) 
confined in a two-dimensional parabolic channel. It was 
shown that the system crystallized in a number of chains 
similar to a colloidal system confined in a narrow chan- 
nel. The authors of Ref. analyzed the structural 
transitions between the ground state configurations and 
showed that such system exhibited a rich phase diagram 
at zero temperature under continuous and discontinuous 
structural transitions. 

The purpose of this study is to achieve a deeper un- 
derstanding of the infiuence of structural properties on 
the motion of paramagnetic colloids in narrow channels 
under non-uniform distributed driving force. In particu- 
lar, we consider: (i) a linear increasing driving force from 
one to the other boundary of the channel, i.e., a linear- 
gradient driving force, and (ii) a parabolic-profile driving 
force, i.e., equal to zero at the boundaries and maximum 
at the center of the channel. Although the experimen- 
tal realization of such profiles is not as obvious as for the 
simplest case of a constant driving force (created, e.g., by 
gravity in an inclined channel), they can serve as a model 
of a very important case of a channel modulated in a non- 
uniform manner, e.g., by a non- uniform distributed dis- 
order ("pinning") or surface roughness. The combination 
of a uniform force applied to colloids with a non-uniform 
modulation in such a channel will result in an effective 
non-uniform driving, i.e., Fdr,nu = Fdr,u - fm,nu- The 
spatial period of the modulation should be chosen, obvi- 
ously, essentially smaller than the colloidal radius. For 
example, if the modulation depth of a channel linearly 
grows from one edge of the channel to another, the effec- 
tive driving force distribution is described by a linearly 
growing function. Moreover, the application of a non- 
uniform driving force to a discrete system of interacting 
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particles confined to a quasi- ID channel is expected to 
result in a rich dynamics interesting from the point of 
view of fundamental research, and it could be useful for 
understanding the dynamics of other interacting systems 
moving in narrow channels. 

The paper is organized as follows. The model is de- 
scribed in Sec. II. In Sec. Ill, we present the results of our 
numerical calculations of the response of the colloidal sys- 
tem to a non- uniform, i.e., linear-gradient and parabolic, 
driving force. In Sec. IV, we study the mobility of stripes 
and dynamical phases of their motion. The conclusions 
are presented in Sec. V. 



II. THE MODEL AND METHOD 

Motion of the system of paramagnetic colloids in a 
two-dimensional (2D) narrow channel under the action 
of a non- uniform distributed driving force Fdrijj) is in- 
vestigated using the Langevin equations of motion in the 
overdamped regime, i.e., the Brownian dynamics (BD) 
method. It is supposed that a magnetic field is ap- 
plied perpendicularly to the channel plane (i.e., perpen- 
dicular to the xy— plane). Following Ref. 0, we choose 
Bo = 1 mT. In this case the magnetic field induces mag- 
netic moments in the colloids Mq = x^o (x = 3 ■ 10^^^ 
0) and, therefore, the interparticle interaction is de- 
scribed by the dipolar repulsive potential V{rij) (for 
fij ^ ^co/, and by hard core interaction for r.y < rcoi)'- 
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where /io is a magnetic constant, fij is the distance be- 
tween ith and jth particles. Following Ref. @, we have 
chosen the colloidal radius rcoi = 2.275 /im. The diame- 
ter of a colloidal particle is denoted as acoi = '2rcoi. 

The use of the BD method assumes that colloids are 
considered in the limit of large viscosity. The system of 
overdamped equations of motion in this case is given by 



' dt 
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where the friction constant t] — 1.2275 • 10~® is Stokes 
friction for a spherical particle of radius Vcoi |3] , is the 
position vector of ith colloid. An external non-uniform 
stationary force FdriVi) is applied along the channel and 
depends only on the transverse coordinate y. The ther- 
mal stochastic term fJ entering Eq. ^ obeys the fol- 
lowing conditions: 



and 
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To find the ground state (or an initial state with the 
lowest free energy close to the ground state) of the sys- 
tem, we first solve Eq. ([2]), for nonzero temperature in 
the absence of driving, thus simulating annealing process 
used for obtaining the ground state of, e.g., colloid s llal 
or vortices in a superconductor (see, e.g., [13, [H, Il9||). 
Then we set temperature equal to zero and solve Eq. ^ 
for colloids driven by the external force FdriUi)- Note 
that inclusion of temperature fluctuation would lead to 
the degradation of the colloidal stripe structure and, as 
a result, to smearing of the velocity profile of different 
colloidal chains. We focus on the features of the veloc- 
ity profiles related to the structure of the chains driven 
by a non-uniform force, and we neglected temperature- 
induced fluctuations in our calculations. In our numerical 
simulations, we use a simulation cell with sizes equal to 
the channel width in the y-direction and typically \AQrcoi 
along the channel, i.e., in the x-direction. We use peri- 
odic boundary condition along the x-direction. Interac- 
tion of colloids with lateral walls is hard- wall. The system 
of overdamped equations ^ is solved using a variable 
time step, calculated from the condition max{dYi\ < S, 
where dr^ is the displacement of ith colloid, and 6 was 
typically chosen rco;/200 or rco;/100, although in most 
cases the value Tcoi/'^O was sufficient. 



III. RESPONSE OF THE SYSTEM TO 
NON-UNIFORM DRIVING FORCE 

A. The ground state configurations 

As was shown in Refs. [1, [ll|, in the absence of ex- 
ternal forces interacting colloids in narrow 2D channels 
are ordered in stripe structures. Therefore, one might 
expect that appliyng an external force along the channel 
will result in an ordered drift of colloidal stripes in the 
channel. (This obviously applies to the case of rather 
weak driving force or its gradient and very narrow chan- 
nels. Strong inhomogeneous driving in a wide enough 
channel could lead to the appearance of instabilities in 
the transverse direction to the applied driving and thus 
to the migration of colloids between stripes.) Thus the 
investigation of the motion of colloids in narrow channels 
can be reduced to the analysis of the drift of colloidal 
stripes. Note that in this case the order of individual 
colloids in a stripe does not change with time, and the 
motion of colloids in s trip es resembles single-file diffusion 
(see, e.g., SIM HI HI Hi). 

For studying the colloidal motion in channels we choose 
such ground state configurations which are characterized 
by strongly pronounced stripe structure and thus by a 
high value of stripe order parameter 
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(4) The order parameter — 1 for particles distributed 
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FIG. 1: (Color online) A typical view of a highly-ordered 
ground-state (^ni = 0.96) colloidal distribution in a channel, 
for the total density p = 0.6(t; 
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B. The integral of motion 

Before discussing results of our numerical calculations, 
we note a useful property that directly follows from sys- 
tem of equations ([2]), namely, the integral of motion: 



dt 
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equidistantly in n; chains across the channel, and ^'„j ^ 
1 for a non-chained case. Here Ly is the channel width, yi 
is the y coordinate of jth colloid and iV is the number of 
particles in the system. Such a choice of the initial con- 
figurations predetermined a well-defined stripe structure 
and allowed us to focus on the study of stripe dynam- 
ics rather than on the dynamics of a low-ordered phase 
(i.e., for '^ni ^ !)■ A typical configuration of the initial 
colloidal distribution is shown in Fig. [1] Our estimate 
of the inter-chain interaction shows that the main con- 
tribution is given only by the adjacent stripes while the 
contribution from more remote stripes is negligible. 

The profile of the applied driving force and the struc- 
tural distribution of colloids in a channel are the two 
most important factors defining the dynamics of colloidal 
motion in a narrow two-dimensional channel in the over- 
damped regime. Their roles essentially distinguish de- 
pending on the value of the applied driving. For an 
overall strong drive the distribution of the driving force 
Fdriy) is dominating. Contrarily, in case of weak drive, 
the motion of colloidal stripes is mainly determined by 
the structural distribution of colloids. In intermediate 
cases the character of colloidal motion results from a com- 
plex competition between these two factors. 

The distribution of the density of colloids in the trans- 
verse direction Pchainiy) is essentially non- uniform, i.e., 
presence of the channel boundaries leads to a relatively 
higher colloidal density on the periphery and lower in the 
central chains ^, i^, ^ . The distribution of colloidal den- 
sity becomes more homogeneous with increasing channel 
width Ly, and in the limit of a wide channel, Ly 3> (TcoI, 
the colloidal structure corresponds to a regular hexagonal 
lattice. 

Another essential feature of the colloidal distribution 
in a channel consists in oscillating dependence of the de- 
fect concentration on the channel width Ly . Thus in- 
crease of the channel width is periodically accompanied 
by a qualitative rearrangement of the colloidal structure, 
namely, by the occurrence of new chains. The values 
of the channel width, for which the colloidal structures 
rearrange, correspond to low-ordered states (^'„, -C 1) 
when the stripe structure is hardly distinguishable. This 
fact is reflected in the oscillation of the stripe order pa- 
rameter VPnj with increasing width of the channel. It is 
worth to note however that the fact that is large does 
not mean that the defect concentration is low. Actually 
it means only that stripes are well pronounced and all 
defects are mainly topological. 



Separating the longitudinal component we obtain, in 
terms of average stripe velocities, 

V^V^n=^Fdr{yn), (7) 

n n 

where n denotes the nth chain, Vxn is the average velocity 
of nth stripe. This expression has a clear meaning: the 
sum of all average chain velocities is defined by the mo- 
mentum of the external force transferred to the system. 
From Eq. ^ it follows that the colloidal system is in the 
static state, i.e., does not move, only in case of: (i) zero 
driving force; (ii) sign-alternating driving force satisfying 
the condition J^n^driVn) = and max \ Fdr{Vn) |< F,,, 
where F^ is some critical force. 

This result is remarkable due to the fact that it differs 
from the case of particles moving under the action of an 
external force in any stationary periodic potential. In the 
latter case, for a commensurate configuration (i.e., when 
the number of colloids per unit length coincides with the 
number of potential minima) , there is always a threshold 
value of the force related to the "static friction". Thus 
the system of particles moves if the driving force exceeds 
the threshold value. In our case, the system of colloids 
moves always when the sum of forces in ((T]) is larger than 
zero. 

Note that the conservation law ([7]) is fulfilled for any 
distribution of the driving force Fdr (y) and it is not valid 
if temperature is non zero. 

C. Stripe velocity versus non-uniform driving force 

In this subsection, we study a response of the colloidal 
system confined in a channel to an external non-uniform 
driving force. Clearly, this case is more general than the 
earlier studied case of a uniform driving created, e.g., by 
the gravity in inclined channels Q. A gradient driving 
force produces a shear stress that allows us to examine 
elastic properties of the colloid "lattice" and to reveal 
the onset of plastic motion. The transition from elastic 
to plastic motion leads to a number of different dynamical 
regimes. 

Without loss of generality, we consider two typical pro- 
files of the non-uniform distributed driving force: (i) a 
linear distributed force given by 

Fdriy) = Fraa. + ' («) 
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FIG. 2: (Color online) The profile of average chain veloc- 
ities {vx) as a function of the transverse coordinate y for a 
parabolic ((a) to (c)) and a linear ((d) to (f)) profiles of the 
applied driving force for various values of a: a = 0.01 to 0.08 
((a), (d)), a = 0.1 to 0.8 ((b), (e)), and a = 1 to 8 ((c), 
(f)). The width of the channel is Ly/ocoi = 4.3 and the total 
colloidal density is p = 0.6 cr^^. The symbols connected by 
solid lines are average velocities of chains while dashed lines 
show velocities of (non-interacting) chains corresponding to 
the distribution of the applied force. 



and (ii) a parabolic distributed driving 




(9) 



where F^ax — mgsin{a) is the maximum value of the 
driving force which we express, for convenience, via the 
projection of the gravity force that acts on a free colloidal 
particle being on an inclined plane at angle a; m is the 
mass of the colloidal particle for the density of colloids 
p,oi = 1600 kg/m^ M- 

The average velocity of different chains function 
of their transverse coordinate y„ is presented in Fig. [2] 
for values of a varying in the range 0.01 < a < 8. The 
profiles of the calculated average stripe velocity Vx(jjn) 
are shown by different symbols (corresponding to differ- 
ent values of driving force expressed in terms of angle a, 
- see the figure captions) connected by solid lines. For 
comparison, the profiles of stripe velocities in the absence 
of interaction between the stripes (i.e., renormalized ap- 
plied external force Fdriy)/^]) are shown by dashed fines, 
for different magnitudes of the driving force. With in- 
creasing a and, accordingly, the total amplitude Fmax of 
the driving force the role of the shape of the distribution 
of applied force on VxiUn) increases, while the influence 



of the structure of colloids and the density of individ- 
ual chains, Pchainiu), on the contrary, decreases. In the 
limit of very large driving force F^r (1 ^ a < 8), the 
distribution of stripe velocity is completely defined by 
the profile of the applied force Fdr{y), as one can ex- 
pect in the plastic limit (see Figs. [21[c),(f), for parabolic 
and linear-gradient driving, correspondingly). On the 
contrary, for small amplitudes Fmax of the driving force 
(0.01 < a ^ 0.05), deviations of Vxn from the profile of 
the driving force Fdr{y) are essential (Fig.[2{a),(d)). This 
limit corresponds to the elastic or quasi-elastic regime. 
In the intermediate region of Fdr{y), 0.05 ^ a < 1, the 
function VxiVn) is essentially infiuenced by both factors, 
Fdriy) and pchamiy)- 

The above behavior is rather general and is typical, in 
principle, for any profile of the external driving. At the 
same time, the shape of the driving-force profile is im- 
portant, especially, in the region of intermediate driving 
forces. In other words, for the same set of parameters 
(i.e., the total density, p, the channel width, Lj,, and a) 
the difference between Fdr and Vx varies for parabolic 
and linear profiles of the external driving force. 

In order to introduce a measure of conformity of the 
profiles of the calculated average velocity Vx{y) and the 
external force Fdr, we define a normalized standard de- 
viation A in the form: 
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Low values of A correspond to small difference between 
Vxniy) and Fdriy)- For a parabolic driving-force profile, 
the gradient of the applied force and thus the shear stress 
between inner chains is much less than between those at 
the periphery, while for a linear profile the gradient is 
constant. 

The results of our calculation of the standard deviation 
A are shown in Fig. [3] versus driving force (a) and density 
(b,c), for linear and parabolic profiles of the driving force. 

Note that for small driving forces, i.e., for elastic mo- 
tion, the deviation of the velocity profile Vxn (y) from that 
of the applied driving Fdriy) is l^ss in case of a linear 
driving-force profile. This means that the elastic defor- 
mation of the colloidal "solid" related to a linear driving 
is stronger than that for a parabolic driving. Contrary, 
for overall large driving (i.e., a), linear-gradient-driven 
stripes follow the profile of the driving force to a lesser 
extent than in case of parabolic driving (see Fig. [3Ja)). 
This means that the average dynamical friction between 
chains is larger in case of a linear-driving applied force. 
This conclusion seems to be quite surprising keeping in 
mind the large gradient of the driving force (experienced 
by peripheral chains) in case of a parabolic driving. As 
we show below, it can be understood in terms of incom- 
mensurability between peripheral chains (i.e., there the 
driving gradient is maximum) related to different densi- 
ties of colloids in different chains. The function A(p) ver- 
sus the overall colloidal density is presented in Figs, ^h) 



5 



<] 0.1 ^ 



(a) 



a 



0.9 
0.8 
0.7 
0.6 
. 0.5 
^ 0.4 
0.3 
0.2 
0.1 
0.0 



0.9 
0.8 
0.7 
0.6 
0.5 
' 0.4 
0.3 
0.2 
0.1 
0.0 



(b) 



•r 



0.3 0.4 0.5 0.6 

-2 



p/cj 



(c) 



0.3 0.4 0.5 0.6 
-2 



p/a 

FIG. 3: (Color online) The normalized standard deviation A 
as a function of a, calculated from data presented in Fig. [2] 
(i.e., for a four-chain structure calculated for p — 0.6 o-'^J. 
Squares (connected by a dot-dashed line) correspond to a 
linear profile of the driving force; circles correspond to a 
parabolic profile (a). The function A vs. the total colloidal 
density p for: a parabolic driving force (b) and a linear- 
gradient driving force (c), for different values of a = 0.1 to 0.8 
(i.e., for three- and four-chain structures as shown in Figs. 4 
and 5). Here we use the same symbols as in Fig. [2] to denote 
velocities of different colloidal stripes. 



and [3jc) , respectively, for parabolic and linear-gradient 
driving, for different a, and Ly/ucoi — 4.3. Note that 
in both cases A(p) has a minimum at about p — O.Su^'^ 
corresponding to a perfect chained structure and incom- 
mensurate number of colloids in different chains. As a 
result, the chains easily slide with respect to each other, 
and the velocity profile reproduces that of the driving 
force (if the driving is strong enough, see Figs, ^ih) and 
(c)). The minimum in A(p) aX p — 0.5(J^^ is followed by 
a sharp increase of A(/9) at about p — O.Gcr^^ explained 
by the appearance of defects in the chained structure. 
The defects lock the motion of different chains, similarly 
to the case of low densities (i.e., p ~ O.Str^'^), and the 
system displays a rigid-body(RB)-like motion. With in- 
creasing p, the maximum in A alternates with a grad- 
ual decrease indicating increasing ordering of colloids in 



stripes. Note that the behavior of A is in general similar 
for the linear and parabolic driving force (i.e., for high 
densities p) , the absolute values and dispersion of A for 
different a are lower in case of the linear-gradient driving, 
in agreement with the result shown in Fig. [3ja). 



D. The role of colloidal density in stripes 

In order to study the influence of two competing fac- 
tors, the structural factor Pchain{y) (i.e., the transverse 
colloidal density) and the distribution of the external 
force Fdr{y), we performed simulations of colloidal mo- 
tion for various values of the total density p in the inter- 
mediate range of values of a, 0.1 < a < 0.8. The average 
stripe velocity function of the transverse coordi- 

nate y at various values of a and different total colloidal 
density p is presented in Figs. |4] and [5l respectively, for 
parabolic and linear-gradient profiles of the external force 
Fdr{y)- The corresponding distributions of the (normal- 
ized) cross-section chain density Pchain{y) are shown in 
the bottom panels of Fig. ID 

As seen from Fig. Ufa), for the density p — O.Slu^'^ and 
a < 0.3 the colloidal structure consists of three stripes 
and moves as a whole (RB motion). If the colloidal dis- 
tribution is symmetric with respect to the channel axis 
y = then the critical value a of the RB-to-plastic mo- 
tion increases up to 0.5. Such an increase of the critical 
value of a is explained by the simultaneous locking of 
motion of the central stripe by two adjacent stripes. The 
transition from the RB mode to the plastic mode occurs 
only for high enough magnitudes of Fmax when colloids 
of the central stripe are capable to overcome the poten- 
tial barriers created by peripheral stripes. However, in 
case of an asymmetric (with respect to the channel axis 
y — Q) distribution Pchain{y), the motion of the central 
stripe is assisted by defects. The potential-energy profile 
asymmetry, together with the transverse degree of free- 
dom, leads to such type of motion of the central stripe 
when it avoids obstacles and moves along the potential- 
energy minimum lines, i.e., along serpentine- like trajec- 
tories. This kind of motion provides a lower, as compared 
to the symmetric case, critical value of a of the transition 
from RB to plastic mode. 

For density p = 0.41, the RB mode is observed only for 
the lowest value a — 0.1 (Fig.lJlJb)). For larger densities 
the RB mode of motion is not observed in the range of 
a =0.1 to 0.8 (Fig.[3]^c)-(f)) but it can be realized for even 



smaller values of a. Note, however, that for p = 0.61cr^; 
(Fig-IHd)) (the least dense state of a four-chain phase) at 
a = 0.1 the motion almost corresponds to the RB mode. 

Comparison of the cases of a parabolic (Fig. 2]) and a 
linear (Fig. [5]) driving forces clearly shows the difference 
in the conditions of the manifestation of the RB mode. 
In particular, for a linear-gradient force Fdr{y) this mode 
is observed for values a < 0.1 and total density p = 
0.31f7^j (i.e., the lowest-p three-chain phase), that means 
a lower threshold magnitude F^ax of the transition from 
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FIG. 4: (Color online) Top panels: The profile of average 
chain velocities {v^) as a function of the transverse coordinate 
y for a parabolic profile of the applied driving force for various 
values of a (a = 0.1 to 0.8) and for the relative width of the 
channel Ly/acoi = 3.7. Bottom panels: The corresponding 
distributions of the chain densities pchainiv) shown for com- 
parison (i.e., normalized on pmax ; where pmax is the maximum 
density of the chain in the colloidal n;-chain configuration). 
The results are shown for the total density p: p — O.SIct^J^^ 
(a), p = 0.41a-^ (a), p = O.Sla"^ (c), p = 0.61aJ (d), 
p = 0.710--^ (e), and p = O.Slcr^^^ (f). As in Fig. [1 symbols 
connected by solid lines show average velocities of chains, and 
dashed lines show velocities of (non-interacting) chains. 



elastic to plastic mode of colloidal motion (cp. Figs.[4l^a) 
andinfa)). At the same time, for the lowest-p four-chain 
phase, this threshold value is higher (a < 0.4) in case 
of a linear driving than for a parabolic, as can be seen 
from Figs.lHd) and [DJd). This means that the RB mode 
manifests itself for: (i) a parabolic driving force Fdr{y) 
in the case of the lowest-p three-chain phase and (ii) a 
linear Fdr{y) in case of the lowest-p four-chain phase. 

The observed behavior is explained as follows. In case 
of a parabolic driving of the lowest-p three-chain struc- 
ture, the driving is applied to the central stripe. The 
motion of this stripe is locked by the potential created 
by two adjacent stripes. The most stable is a symmet- 
ric configuration, i.e., when two peripheral stripes are 
mirror-symmetric with respect to the axis of the channel 
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FIG. 5: (Color online) The profile of average chain veloc- 
ities {vx) as a function of the transverse coordinate y for 
a liner-gradient profile of the applied driving force for vari- 
ous values of a (a =0.1 to 0.8) and for relative width of the 
channel Ly/acoi = 3.7. The results are shown for the total 
density p: p = 0.31(7"^ (a), p = 0.41cr-^ (a), p = 0.51(7"^ 
(c), p = 0.61a-^ (d), p = 0.71a-^ (e), and p = 0.81a J (f). 
As in Fig. [2] symbols connected by solid lines show average 
velocities of chains, and dashed lines show velocities of (non- 
interacting) chains. 



?/ = 0, and thus the potential created by them is double 
the potential created by a single stripe. Thus, the sym- 
metric configuration favors the RB mode. For asymmet- 
ric configurations, the RB-to-plastic mode threshold is 
lower. In this case the motion of the central stripe (i.e., 
for driving higher than a critical value) is serpentine- 
like, i.e., the central stripe deforms in the y-direction 
following the minimum potential-energy path. A typical 
dispersion of the trajectory of the central stripe in the 
transverse direction is Ay « 0.1 (while Ay « 10~* to 
10"'^ for the symmetric case). As a result, the lowest-p 
three-chain colloidal structure is characterized by a rela- 
tively high threshold of the RB-to-plastic transition, un- 
der a parabolic-profile driving (see Fig.^a)). In case of a 
linear-gradient driving of the lowest-p three-chain phase 
(Fig-EJa)), the maximum driving force is applied to one 
of the peripheral chains. In the RB mode, its motion 
is locked only by the adjacent central chain. When un- 
locked, the peripheral chain also provides an additional 
driving applied to the central chain. The central chain 
adjusts itself to follow the minimum potential-energy 
path between the two peripheral chains. This motion 
is characterized by a rather large transverse dispersion. 
Ay « 0.2 to 0.3, that facilitates an easy sliding of the 
three chains with respect to each other. 
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FIG. 6: (Color online) Average stripe velocities {vx) as a 
function of the total colloidal density p for a linear ((a), (b)) 
and a parabolic ((c), (d)) profiles of the applied driving force, 
for the channel width Ly/ocoi = 5, and for different driving: 
a = 0.3 ((a), (c)) and a = 0.7 ((b), (d)). Different symbols 
connected by solid lines show (vx) of different chains marked 
for clarity by numbers f, 2, 3, 4, and 5 (the numbers are 
shown by the same color as the corresponding curves). 



The situation is quite different in case of the lowest-p 
four-layer phase shown in Figs. Hljd) and O^d). In the 
ground-state low-p four-chain configurations, two central 
chains have lower density than the peripheral chains (see 
bottom panels in Figs. El^d), (e), and (f)). Thus, the 
potential profile created by the central chains is deeper, 
and they survive a rather strong shear stress (similar to 
the case of low-density ID chains in a harmonic potential 
in the Frenkel-Kontorova model [25*]). Such shear stress 
is negligible for a parabolic driving (it appears as a re- 
sult of transverse instabilities). Thus the central part 
of the four-chain structure remain in the RB mode for 
very large values of driving, and it also does in case of 
a linear-gradient driving (see Figs. [Jd), (e), and (f)). 
However, the higher-p peripheral chains produce shal- 
lower potential-energy profiles, and thus the friction be- 
tween them and the central chains is weak [2^ . It is clear 
than that in case of a parabolic driving the peripheral 
stripes can easily slide with respect to the central stripes, 
and thus the RB-to-plastic mode threshold is rather low 
(Fig. Ud)). In contrast, a relatively weak (i.e., for the 
same Fmax) gradient of the driving force in case of a 
linear-gradient driving provides a high threshold of the 
RB-to-plastic mode transition (Fig.lSJd)). 

As we discussed above, the colloidal distribution in a 
channel depends on the width of the channel and the to- 
tal colloidal density. For certain "matching" sets of these 
parameters, colloids form well-defined stripe structures 
characterized by a maximum value of the order param- 
eter = 1, while small values of the order parameter 
1 correspond to non-chained structure. It is in- 
teresting then to follow the evolution of stripe velocities 



when changing the total colloidal density in the channel 
[26| . The results of calculations are presented in Fig. [5] 
for a linear-gradient ((a), (b)) and a parabolic ((c), (d)) 
profile of the driving force and different magnitudes of 
driving. Let us consider first the case of a linear driving. 
At low density, p — O.Scr^j, the colloids form four stripes 
(marked by numbers 1, 2, 3, and 4 of the same color as 
the corresponding lines/symbols), and for all values of a 
(i.e., a = 0.3 (a) and 0.7 (b)) shown in Fig.[6l two central 
chains (2 and 3) turn out to be locked while the periph- 
eral chains (1 and 4) (which have a higher stripe density 
Pchain and are incommensurate with the central chains) 
slide with respect to the central chains. With increas- 
ing p, e.g., p = 0.4cr^^ for a — 0.7 (see Fig. [SJb)) the 
central chains 2 and 3 unlock. This is explained by de- 
creasing of the inter-colloid distance in stripes and thus 
by shallowing of the potential profile created by them. 
Note that for a weaker driving a — 0.3 (a), the central 
stripes remain locked at p = ^■^'^coi they unlock at 
a higher value of p: p = O.Scr^'^. For even higher den- 
sity, p = 0.6cr^j, the inter-shell defects lock the motion 
of the central chains 2 and 3 (and partially of all the 
chains), and the whole stripe structure collapses. If we 
still increase the colloid density, a new stripe (5) origi- 
nates from the "disordered" phase, and for larger p (i.e., 
p — 0-8a^^i) the motion of all five stripes unlock, and 
they move with individual velocities. Note that the tran- 
sition of the system from the state with ni stripes to the 
state with ni + 1 stripes occurs always through the dis- 
ordered phase (i.e., when all or some of the stripes lose 
their identity, and colloids of these "stripes" move with 
the same velocity as, e.g., stripes 2, 3, and 4 in Fig. [51^ a) 
and (b) at p = O.Gcr^j). Thus the disordered phase of 
the motion serves as a "bifurcation point" that gives rise 
to a new state with ni + 1 velocity branches. The evolu- 
tion of stripe velocities for a parabolic driving occurs in a 
similar way (Fig. [S{c) and (d)), although the degenerate 
velocities of the central stripes, 2 and 3 (and those of the 
peripheral stripes, 1 and 4, which are separated from the 
velocities of the central stripes by a wide gap), unlock 
only for rather high density. 

It is also interesting to note an oscillating behavior of 
the velocity of motion of the peripheral stripes driven 
only by the interaction with adjacent stripes. The os- 
cillation in velocity refiects the oscillation of the friction 
force between the stripes with increasing total colloidal 
density. 



IV. MOBILITY OF STRIPES AND 
DYNAMICAL PHASES OF THEIR MOTION 

As shown above, the response of the system of col- 
loidal stripes (i.e., the profile of individual velocities of 
different stripes) to the external non-uniform driving is, 
in general, nonlinear: the velocity profile of the stripes 
is not scalable with the driving (i.e., for all values of the 
driving, although there are regions of scalability which 



will be discussed below). 

In order to systematically examine the response of the 
colloidal system to an external driving, we analyze the 
mobility of different chains under the action of a delta- 
like driving force applied to one of the chains: 



(11) 



where S{y - yk,i) = 1, if y = yuA, and 5{y - yk,i) = 0, 
if 2/ 7^ Wfejij yk,i is the coordinate of ith particle of fcth 
stripe [27| and Fmax is the driving force as defined in 
Eq. ([9]). We simulated the colloidal motion in a chan- 
nel of width Ly/acoi = 4 for various colloidal densities p. 
For such a channel and considered densities p, the ground 
state of the system corresponds to a three-stripe struc- 
ture characterized by a high value of the order parameter 
^'3 « 1. Fig. [7] shows the results of calculation for the 
mobilities versus driving force Fdrk for different stripes 
Bn,k = Vxri/Fdrk, whcrc Vxn is the average velocity of 
nth stripe, and Fdrk is the force applied to kih. stripe. 
(Here we show only the case when driving is applied to 
one peripheral stripe. Applying driving to other stripes 
results in a similar behavior.) Note that the mobility is 
defined by B^^k = v^k/Fdrk, while B„^fc for k ^ n is a. 
mobile response of the stripes. 

As seen from Fig. [3 the mobility as a function of the 
applied force exhibits a variety of modes of motion. First, 
we stress the most general features of the behavior. A 
stripe mobility Bk^k tends to unity in the limit of large ap- 
plied force Fdrk, and mobile responses Bn.k, correspond- 
ingly, vanish, as seen from Fig. [71 This limit corresponds 
to the plastic mode of motion. Another general conclu- 
sion refers to the RB-to-plastic transition: this transition 
shifts towards lower drivings with increasing the total 
colloidal density (cp. Figs. [TJa) and (b)). As we dis- 
cussed above, a growth of the colloidal density results in 
a hardening of the rigidity of the elastic stripes [l^ and 
a simultaneous weakening of the inter-stripe friction. 

However, apart from these general features, the mo- 
bilities display other important peculiarities. To discuss 
these, let us define dynamical phases corresponding to 
different regimes of colloidal motion in channels. We 
denote the first regime as "RB" ("rigid body", or elas- 
tic motion). This regime is most pronounced for low 
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(Fig. EKa)), where it ex- 

-2 



colloidal density, p 

tends to Fdri « 1.05. For p = O-bba^^i (Fig. EKb)) the 
threshold value reduced by a factor of five, Fdri ~ 0.2. 
For higher densities, the RB-to-plastic mode transition 
shifts towards very low values of Fdri- The RB regime 
is followed by a broad transition region ( "Tr" ) character- 
ized by a monotonic increase of the mobility Bk,k (and 
monotonic decrease of the mobile responses Bn,k)- While 
this growth is close to linear in Fig. EKa) , for higher den- 
sity it displays oscillations in Fig. [7{b). Remarkably, for 
higher colloidal densities these oscillations turn to regions 
of decreasing colloid mobility Bf^ i^ versus driving force. 
Note that this regime (denoted as "S") is similar to the 
dynamically-ordered phase of vortex motion in supercon- 
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FIG. 7: (Color online) The mobility Bk.k and mobile response 
Bn,k{n 7^ k) of different stripes as a function of driving force 
applied to fcth stripe (fc = 1) at various values of the total 
density: p = 0.45aJ„^ (a), p = CSSa'^ (b), p = CGSa'^ 
(c), and p = 0.75a^Q( (d). The mobility (mobile response) 
of different stripes are shown by different symbols; _Bi,i (□), 
Bi.2 (O), and Bi,3 (A). 



ductors with arrays of regular pinning sites |28l . |29[ . This 
ordered phase follows a disordered phase with a higher 
density of average vortex flow and thus leads to the ap- 
pearance of a negative differential resistivity (NDR) part 
of N-type in the VTcurve [29*1. In our case of colloidal 
motion in a narrow channel, we observe a similar effect 
characterized by a slowing down of the motion of the 
stripe driven by an external force (with a simultaneous 
increase of the velocity of the adjacent stripes) when in- 
creasing the applied driving force. In this region, the 
system of colloids displays a partial "reentrant" behav- 
ior (cp. Refs. HO, [Ml) when, being melted by increasing 
shear stress, it starts evolving towards a solid phase with 
further increasing driving. Note that this dynamically- 
induced "solidification" has common features with the re- 
cently discovered transition "freezing by heating" [s^, HO] 
where by increasing temperature it leads to the crystal- 
lization of moving repulsive interacting particles which 
are in a molten state. 

The mechanism of the observed colloid solidification 
can be understood as follows. As we discussed above, 
the longitudinal motion of colloidal stripes is accompa- 
nied by transverse oscillations of the stripes' trajectories 
(i.e., serpentine-like motion) related to the asymmetry 
of the potential-energy profiles created by the adjacent 
stripes. For low drivings (but larger than the critical 
driving force of the RB-to-plastic transition), the stripe 
driven by an external force moves along some serpentine- 
like trajectory. When moving, the stripe itself deforms 
adjusting to the potential-energy landscape. In its turn, 
it elastically deforms other stripes due to the interaction 
with colloids in the adjacent stripes. At low velocities, 
colloids in adjacent stripes relax to the initial state thus 
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providing low friction between the stripes. Increasing 
the velocity of motion of the driven stripe (at high col- 
loid density) leads to increasing rate of the inter-stripe 
collisions and to the development of instabilities in the 
transverse direction. This results in an increase of the 
dynamical friction between stripes which explains the ob- 
served decrease of the mobility B^^k (and increase of the 
mobile response Bn^k), i-e., the appearance of "S" phase. 
Very large driving, however, straightens the trajectories 
of motion of the driven stripe. In the limit of F^r ~^ oOj 
the driven stripe moves along unrelaxed in the y-direction 
straight trajectory, and the mobility Bk^k — 1 (while the 
mobile response Bn,k 0), as shown in Figs.[7Kc), (d). 

V. CONCLUSIONS 

We have studied the dynamics of colloids driven by 
an external non-uniform force in a narrow channel. We 
focused on colloidal densities which provide well-defined 
colloidal stripe structures. In the limit of a very low 
driving force the system moves as a rigid body (elastic 
motion) independent of the profile of the driving force. 
In the opposite limit of a very strong driving force, the 
profile of the average velocities of colloidal stripes fol- 
low the profile of the applied driving (plastic regime). 
While these limiting-case results are easily understand- 
able, the dynamics of colloidal stripes in the general case 
is rather complex. It is governed by the interplay of sev- 
eral factors, i.e.: (i) the magnitude and profile of driv- 
ing force; (ii) the total density of the colloidal system 
and the distribution of the colloidal density in stripes, 
and, as a consequence, (iii) commensurability effects be- 
tween the adjacent chains. We have shown that depend- 



ing on the density and number of stripes, the transition 
from elastic to plastic motion occurs at different values of 
the driving force. For example, the lowest-density three- 
stripe colloidal configuration is shown to be more ro- 
bust with respect to a parabolic driving than to a linear- 
gradient driving, while for the lowest-density four-stripe 
configuration the situation is opposite. This result is ex- 
plained by the colloidal density distribution over the cen- 
tral (lower density) and the peripheral (higher density) 
stripes and commensurability effects. We have analyzed 
the mobility of colloidal stripes and have identified the 
dynamical phases of their motion. In particular, it has 
been shown that the transition from elastic (rigid-body- 
motion phase) to plastic mode, depending on the density 
of colloids, could be either monotonic (i.e., characterized 
by a gradual increase of the mobility), or it could contain 
an NDR-type part (i.e., characterized by a drop of mo- 
bility versus driving force). This unusual "solidification" 
is related to the dynamically-induced increase of the fric- 
tion in the colloidal system and is similar to the recently 
discovered "freezing by heating" transition. The results 
of our study, with corresponding changes, can also be ap- 
plied to other systems of interacting particles driven by 
a non- uniform force in narrow channels, e.g., in physics 
or biology. 
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